Computational Methods for Detecting Large-Scale Chromosome Rearrangements in SNP Data

Large-scale chromosome rearrangements such as copy number variants (CNVs) and inversions encompass a considerable proportion of the genetic variation between human individuals. In a number of cases, they have been closely linked with various inheritable diseases. Single-nucleotide polymorphisms (SNPs) are another large part of the genetic variance between individuals. They are also typically abundant and their measuring is straightforward and cheap. This thesis presents computational means of using SNPs to detect the presence of inversions and deletions, a particular variety of CNVs. Technically, the inversion-detection algorithm detects the suppressed recombination rate between inverted and non-inverted haplotype populations whereas the deletion-detection algorithm uses the EM-algorithm to estimate the haplotype frequencies of a window with and without a deletion haplotype. As a contribution to population biology, a coalescent simulator for simulating inversion polymorphisms has been developed. Coalescent simulation is a backward-in-time method of modelling population ancestry. Technically, the simulator also models multiple crossovers by using the Counting model as the chiasma interference model. Finally, this thesis includes an experimental section. The aforementioned methods were tested on synthetic data to evaluate their power and specificity. They were also applied to the HapMap Phase II and Phase III data sets, yielding a number of candidates for previously unknown inversions, deletions and also correctly detecting known such rearrangements.

LKB1 tumor suppression in Peutz-Jeghers Syndrome

Kasvainten, ajatellaan syntyvän yksittäisen solun perimän mutaatioista, jonka seurauksena tuon solun kasvu häiriintyy. Ruoansulatuskanavan polyyppien syntyä käytetään usein mallina siitä, miten nämä epiteelisoluun kerääntyvät mutaatiot aiheuttavat asteittain pahenevan kasvuhäiriön. Peutz–Jeghersin oireyhtymä (PJS) on perinnöllinen polypoosisyndrooma, jossa oireita aiheuttavat erityisesti maha-suolikanavan hamartomatoottiset polyypit. Noin puolella PJS potilaista havaitaan mutaatioita LKB1 kasvunrajoite geenissä. Hiirille joilta toinen Lkb1 alleeli on poistettu (Lkb1+/-) kehittyy PJS-tyypin maha-suolikanavan polyyppeja, joissa on epiteelin liikakasvun lisäksi merkittävä sileälihaskomponentti, aivan kuten PJS polyypeissa. Kuten myös muissa ruoansulatuskanavan polypooseissa, sekä PJS että hiirten polyypeissa Cyclo-oxygenaasi-2:n (COX-2) määrä on usein kohonnut. PJS-polyyppien kehittymisen molekulaarinen mekanismi on kuitenkin selvittämättä. Koska vain osa PJS potilaista kantaa LKB1 mutaatioita, mutaatiot jossakin toisessa lokuksessa saattaisivat selittää osan PJS tapauksista. Jotta PJS:n geneettinen tausta selviäisi, seulottiin kolmen LKB1:n kanssa interaktoivan proteiinin (BRG1, STRADα ja MO25α) geenit PJS potilaista joilla ei ole havaittu LKB1 mutaatioita. Yhdessäkään tutkituista geeneistä ei havaittu tautia aiheuttavia mutaatioita. Näiden kolmen geenin pois sulkeminen, ja uusien menetelmien ansiosta kasvanut havaittujen Lkb1 mutaatioden määrä viittaavat LKB1:n olevan useimpien PJS tapausten taustalla. COX-2:n estäjien käyttö on tehokkaasti vähentänyt polyyppien määrää familiaarisessa adenomatoottisessa polypoosissa. Tästä johtuen COX-2:n eston tehokkuutta tutkittiin PJS polypoosissa. PJS-tyypin polypoosin havaittin pienenevän merkittävästi Lkb1+/- hiirissä, joilta oli lisäksi poistettu toinen tai molemmat COX-2:n alleeleista. Lisäksi farmakologinen COX-2:n esto Celecoxib:lla vähensi polypoosia tehokkaasti. Näin ollen COX-2:n eston tehokkuutta tutkittiin seuraavaksi PJS potilaissa. Kuuden kuukauden Celecoxib hoidon jälkeen polypoosin havaittiin vähentyneen merkittävästi osalla potilaista (2/6). Nämä tulokset osoittavat COX-2:n roolin PJS-polyyppien kehityksessä, ja viittaavat COX-2:n eston vähentävän polypoosia. Kasvunrajoitegeenin klassisen määritelmän mukaan kasvaimen kehitys vaatii perinnöllisen mutaation lisäksi geenin toisenkin alleelin mutaation, mutta PJS-polyyppien häiriintyneestä epiteelistä ei kuitenkaan systemaattisesti löydy toista LKB1:n mutaatiota. Havainto johti tutkimukseen, jossa selvitettiin voisiko LKB1:n kasvun rajoitus välittyäkin epäsuorasti tukikudokseksi ajatelluista sileälihassoluista. Tätä tutkittiin kehittämällä poistogeeninen hiirimalli jossa Lkb1 on mutatoitunut vain sileälihassoluissa. Näille hiirille kehittyi polyyppeja, jotka ovat kaikin tavoin PJS-polyyppien kaltaisia. Lkb1:n menettäneiden solujen havaittiin tuottavan vähemmän transformoivaa kasvutekijä beetaa (TGFß), joka aiheutti solujen välisen viestinnän heikentymisen ja mahdollisesti viereisten epiteelisolujen liikakasvun. Vastaava häiriö havaittiin myös PJS-potilaiden polyypeissa, mikä viittaa siihen, että potilaillakin sileälihassolujen häiriö on polyyppien taustalla. Havainto suuntaa täten hoitokohteiden etsintää ja osoittaa että LKB1 toimii kasvunrajoittajana epätyypillisellä tavalla pitäen naapurisolujen kasvun kurissa.

Molecular mechanisms of cancer predisposition in HNPCC/Lynch syndrome

Hereditary non-polyposis colorectal carcinoma (HNPCC; Lynch syndrome) is among the most common hereditary cancers in man and a model of cancers arising through deficient DNA mismatch repair (MMR). It is inherited in a dominant manner with predisposing germline mutations in the MMR genes, mainly MLH1, MSH2, MSH6 and PMS2. Both copies of the MMR gene need to be inactivated for cancer development. Since Lynch syndrome family members are born with one defective copy of one of the MMR genes in their germline, they only need to acquire a so called second hit to inactivate the MMR gene. Hence, they usually develop cancer at an early age. MMR gene inactivation leads to accumulation of mutations particularly in short repeat tracts, known as microsatellites, causing microsatellite instability (MSI). MSI is the hallmark of Lynch syndrome tumors, but is present in approximately 15% of sporadic tumors as well. There are several possible mechanisms of somatic inactivation (i.e. the second hit ) of MMR genes, for instance deletion of the wild-type copy, leading to loss of heterozygosity (LOH), methylation of promoter regions necessary for gene transcription, or mitotic recombination or gene conversion. In the Lynch syndrome tumors carrying germline mutations in the MMR gene, LOH was found to be the most frequent mechanism of somatic inactivation in the present study. We also studied MLH1/MSH2 deletion carriers and found that somatic mutations identical to the ones in the germline occurred frequently in colorectal cancers and were also present in extracolonic Lynch syndrome-associated tumors. Chromosome-specific marker analysis implied that gene conversion, rather than mitotic recombination or deletion of the respective gene locus accounted for wild-type inactivation. Lynch syndrome patients are predisposed to certain types of cancers, the most common ones being colorectal, endometrial and gastric cancer. Gastric cancer and uroepithelial tumors of bladder and ureter were observed to be true Lynch syndrome tumors with MMR deficiency as the driving force of tumorigenesis. Brain tumors and kidney carcinoma, on the other hand, were mostly MSS, implying the possibility of alternative routes of tumor development. These results present possible implications in clinical cancer surveillance. In about one-third of families suspected of Lynch syndrome, mutations in MMR genes are not found, and we therefore looked for alternative mechanisms of predisposition. According to our results, large genomic deletions, mainly in MSH2, and germline epimutations in MLH1, together explain a significant fraction of point mutation-negative families suspected of Lynch syndrome and are associated with characteristic clinical and family features. Our findings have important implications in the diagnosis and management of Lynch syndrome families.

Genotyping for genetic association studies: Methods and applications

In this thesis, two separate single nucleotide polymorphism (SNP) genotyping techniques were set up at the Finnish Genome Center, pooled genotyping was evaluated as a screening method for large-scale association studies, and finally, the former approaches were used to identify genetic factors predisposing to two distinct complex diseases by utilizing large epidemiological cohorts and also taking environmental factors into account. The first genotyping platform was based on traditional but improved restriction-fragment-length-polymorphism (RFLP) utilizing 384-microtiter well plates, multiplexing, small reaction volumes (5 µl), and automated genotype calling. We participated in the development of the second genotyping method, based on single nucleotide primer extension (SNuPeTM by Amersham Biosciences), by carrying out the alpha- and beta tests for the chemistry and the allele-calling software. Both techniques proved to be accurate, reliable, and suitable for projects with thousands of samples and tens of markers. Pooled genotyping (genotyping of pooled instead of individual DNA samples) was evaluated with Sequenom s MassArray MALDI-TOF, in addition to SNuPeTM and PCR-RFLP techniques. We used MassArray mainly as a point of comparison, because it is known to be well suited for pooled genotyping. All three methods were shown to be accurate, the standard deviations between measurements being 0.017 for the MassArray, 0.022 for the PCR-RFLP, and 0.026 for the SNuPeTM. The largest source of error in the process of pooled genotyping was shown to be the volumetric error, i.e., the preparation of pools. We also demonstrated that it would have been possible to narrow down the genetic locus underlying congenital chloride diarrhea (CLD), an autosomal recessive disorder, by using the pooling technique instead of genotyping individual samples. Although the approach seems to be well suited for traditional case-control studies, it is difficult to apply if any kind of stratification based on environmental factors is needed. Therefore we chose to continue with individual genotyping in the following association studies. Samples in the two separate large epidemiological cohorts were genotyped with the PCR-RFLP and SNuPeTM techniques. The first of these association studies concerned various pregnancy complications among 100,000 consecutive pregnancies in Finland, of which we genotyped 2292 patients and controls, in addition to a population sample of 644 blood donors, with 7 polymorphisms in the potentially thrombotic genes. In this thesis, the analysis of a sub-study of pregnancy-related venous thromboses was included. We showed that the impact of factor V Leiden polymorphism on pregnancy-related venous thrombosis, but not the other tested polymorphisms, was fairly large (odds ratio 11.6; 95% CI 3.6-33.6), and increased multiplicatively when combined with other risk factors such as obesity or advanced age. Owing to our study design, we were also able to estimate the risks at the population level. The second epidemiological cohort was the Helsinki Birth Cohort of men and women who were born during 1924-1933 in Helsinki. The aim was to identify genetic factors that might modify the well known link between small birth size and adult metabolic diseases, such as type 2 diabetes and impaired glucose tolerance. Among ~500 individuals with detailed birth measurements and current metabolic profile, we found that an insertion/deletion polymorphism of the angiotensin converting enzyme (ACE) gene was associated with the duration of gestation, and weight and length at birth. Interestingly, the ACE insertion allele was also associated with higher indices of insulin secretion (p=0.0004) in adult life, but only among individuals who were born small (those among the lowest third of birth weight). Likewise, low birth weight was associated with higher indices of insulin secretion (p=0.003), but only among carriers of the ACE insertion allele. The association with birth measurements was also found with a common haplotype of the glucocorticoid receptor (GR) gene. Furthermore, the association between short length at birth and adult impaired glucose tolerance was confined to carriers of this haplotype (p=0.007). These associations exemplify the interaction between environmental factors and genotype, which, possibly due to altered gene expression, predisposes to complex metabolic diseases. Indeed, we showed that the common GR gene haplotype associated with reduced mRNA expression in thymus of three individuals (p=0.0002).

The Myopathic Protein Myotilin in Developing Mouse and in Muscle Function

Skeletal muscle cells are highly specialised in order to accomplish their function. During development, the fusion of hundreds of immature myoblasts creates large syncytial myofibres with a highly ordered cytoplasm filled with packed myofibrils. The assembly and organisation of contractile myofibrils must be tightly controlled. Indeed, the number of proteins involved in sarcomere building is impressive, and the role of many of them has only recently begun to be elucidated. Myotilin was originally identified as a high affinity a-actinin binding protein in yeast twohybrid screen. It was then found to interact also with filamin C, actin, ZASP and FATZ-1. Human myotilin is mainly expressed in striated muscle and induces efficient actin bundling in vitro and in cells. Moreover, mutations in myotilin cause different forms of muscle disease, now collectively known as myotilinopathies. In this thesis, consisting of three publications, the work on the mouse orthologue is presented. First, the cloning and molecular characterisation of the mouse myotilin gene showed that human and mouse myotilin share high sequence homology and a similar expression pattern and gene regulation. Functional analysis of the mouse promoter revealed the myogenic factor-binding elements that are required for myotilin gene transcription. Secondly, expression of myotilin was studied during mouse embryogenesis. Surprisingly, myotilin was expressed in a wide array of tissues at some stages of development; its expression pattern became more restricted at perinatal stages and in adult life. Immunostaining of human embryos confirmed broader myotilin expression compared to the sarcomeric marker titin. Finally, in the third article, targeted deletion of myotilin gene in mice revealed that it is not essential for muscle development and function. These data altogether indicate that the mouse can be used as a model for human myotilinopathy and that loss of myotilin does not alter significantly muscle structure and function. Therefore, disease-associated mutant myotilin may act as a dominant myopathic factor.

Genetic basis of hereditary colorectal cancers: Hereditary nonpolyposis colorectal cancer and Familial adenomatous polyposis

Hereditary nonpolyposis colorectal cancer (HNPCC) and familial adenomatous polyposis (FAP) are characterized by a high risk and early onset of colorectal cancer (CRC). HNPCC is due to a germline mutation in one of the following MMR genes: MLH1, MSH2, MSH6 and PMS2. A majority of FAP and attenuated FAP (AFAP) cases are due to germline mutations of APC, causing the development of multiple colorectal polyps. To date, over 450 MMR gene mutations and over 800 APC mutations have been identified. Most of these mutations lead to a truncated protein, easily detected by conventional mutation detection methods. However, in about 30% of HNPCC and FAP, and about 90% of AFAP families, mutations remain unknown. We aimed to clarify the genetic basis and genotype-phenotype correlation of mutation negative HNPCC and FAP/AFAP families by advanced mutation detection methods designed to detect large genomic rearrangements, mRNA and protein expression alterations, promoter mutations, phenotype linked haplotypes, and tumoral loss of heterozygosity. We also aimed to estimate the frequency of HNPCC in Uruguayan CRC patients. Our expression based analysis of mutation negative HNPCC divided these families into two categories: 1) 42% of families linked to the MMR genes with a phenotype resembling that of mutation positive, and 2) 58% of families likely to be associated with other susceptibility genes. Unbalanced mRNA expression of MLH1 was observed in two families. Further studies revealed that a MLH1 nonsense mutation, R100X was associated with aberrant splicing of exons not related to the mutation and an MLH1 deletion (AGAA) at nucleotide 210 was associated with multiple exon skipping, without an overall increase in the frequency of splice events. APC mutation negative FAP/AFAP families were divided into four groups according to the genetic basis of their predisposition. Four (14%) families displayed a constitutional deletion of APC with profuse polyposis, early age of onset and frequent extracolonic manifestations. Aberrant mRNA expression of one allele was observed in seven (24%) families with later onset and less frequent extracolonic manifestations. In 15 (52%) families the involvement of APC could neither be confirmed nor excluded. In three (10%) of the families a germline mutation was detected in genes other than APC: AXIN2 in one family, and MYH in two families. The families with undefined genetic basis and especially those with AXIN2 or MYH mutations frequently displayed AFAP or atypical polyposis. Of the Uruguayan CRC patients, 2.6% (12/461) fulfilled the diagnostic criteria for HNPCC and 5.6% (26/461) were associated with increased risk of cancer. Unexpectedly low frequency of molecularly defined HNPCC cases may suggest a different genetic profile in the Uruguayan population and the involvement of novel susceptibility genes. Accurate genetic and clinical characterization of families with hereditary colorectal cancers, and the definition of the genetic basis of "mutation negative" families in particular, facilitate proper clinical management of such families.

A capability-based view of organisational renewal : Combining opportunity- and advantage-seeking growth in large, established European and North American wood-industry companies

The purpose of this study was to extend understanding of how large firms pursuing sustained and profitable growth manage organisational renewal. A multiple-case study was conducted in 27 North American and European wood-industry companies, of which 11 were chosen for closer study. The study combined the organisational-capabilities approach to strategic management with corporate-entrepreneurship thinking. It charted the further development of an identification and classification system for capabilities comprising three dimensions: (i) the dynamism between firm-specific and industry-significant capabilities, (ii) hierarchies of capabilities and capability portfolios, and (iii) their internal structure. Capability building was analysed in the context of the organisational design, the technological systems and the type of resource-bundling process (creating new vs. entrenching existing capabilities). The thesis describes the current capability portfolios and the organisational changes in the case companies. It also clarifies the mechanisms through which companies can influence the balance between knowledge search and the efficiency of knowledge transfer and integration in their daily business activities, and consequently the diversity of their capability portfolio and the breadth and novelty of their product/service range. The largest wood-industry companies of today must develop a seemingly dual strategic focus: they have to combine leading-edge, innovative solutions with cost-efficient, large-scale production. The use of modern technology in production was no longer a primary source of competitiveness in the case companies, but rather belonged to the portfolio of basic capabilities. Knowledge and information management had become an industry imperative, on a par with cost effectiveness. Yet, during the period of this research, the case companies were better in supporting growth in volume of the existing activity than growth through new economic activities. Customer-driven, incremental innovation was preferred over firm-driven innovation through experimentation. The three main constraints on organisational renewal were the lack of slack resources, the aim for lean, centralised designs, and the inward-bound communication climate.

Eurasian Neogene large herbivorous mammals and climate

The development and changes in the distribution of herbivorous mammal communities during the Neogene is complex. The Eurasian scale environmental patterns reflect the large scale geographical and climatic patterns. The reorganization of these affect the biome distribution throughout the continent. The distribution of mammal taxa was closely associated with the distribution of biomes. In Eurasia the Neogene development of environments was twofold. The early and middle Miocene that seemed to have been advantageous for mammals was followed by drying of environments during the late Neogene. The mid-latitude drying was the main trend, and it is the combined result of the retreat of Paratethys, the uplift of Tibetan Plateau and changes in the ocean currents and temperatures. The common mammals were "driving" the evolution of mammalian communities. During the late Miocene we see the drying affecting more and more regions, and we see changes in the composition of mammalian communities.

Forest inventory-based large-scale forest biomass and carbon budget assessment: new enhanced methods and use of remote sensing for verification

In recent years, concern has arisen over the effects of increasing carbon dioxide (CO2) in the earth's atmosphere due to the burning of fossil fuels. One way to mitigate increase in atmospheric CO2 concentration and climate change is carbon sequestration to forest vegeta-tion through photosynthesis. Comparable regional scale estimates for the carbon balance of forests are therefore needed for scientific and political purposes. The aim of the present dissertation was to improve methods for quantifying and verifying inventory-based carbon pool estimates of the boreal forests in the mineral soils. Ongoing forest inventories provide a data based on statistically sounded sampling for estimating the level of carbon stocks and stock changes, but improved modelling tools and comparison of methods are still needed. In this dissertation, the entire inventory-based large-scale forest carbon stock assessment method was presented together with some separate methods for enhancing and comparing it. The enhancement methods presented here include ways to quantify the biomass of understorey vegetation as well as to estimate the litter production of needles and branches. In addition, the optical remote sensing method illustrated in this dis-sertation can be used to compare with independent data. The forest inventory-based large-scale carbon stock assessment method demonstrated here provided reliable carbon estimates when compared with independent data. Future ac-tivity to improve the accuracy of this method could consist of reducing the uncertainties regarding belowground biomass and litter production as well as the soil compartment. The methods developed will serve the needs for UNFCCC reporting and the reporting under the Kyoto Protocol. This method is principally intended for analysts or planners interested in quantifying carbon over extensive forest areas.

Ecological processes and large-scale climate relationships in northern coniferous forests

Ilmasto vaikuttaa ekologisiin prosesseihin eri tasoilla. Suuren mittakaavan ilmastoprosessit, yhdessä ilmakehän ja valtamerien kanssa, säätelevät paikallisia sääilmiöitä suurilla alueilla (mantereista pallopuoliskoihin). Tämä väistöskirja pyrkii selittämään kuinka suuren mittakaavan ilmasto on vaikuttanut tiettyihin ekologisiin prosesseihin pohjoisella havumetsäalueella. Valitut prosessit olivat puiden vuosilustojen kasvu, metsäpalojen esiintyminen ja vuoristomäntykovakuoriaisen aiheuttamat puukuolemat. Suuren mittakaavan ilmaston löydettiin vaikuttaneen näiden prosessien esiintymistiheyteen, kestoon ja levinneisyyteen keskeisten sään muuttujien välityksellä hyvin laajoilla alueilla. Tutkituilla prosesseilla oli vahva yhteys laajan mittakaavan ilmastoon. Yhteys on kuitenkin ollut hyvin dynaaminen ja muuttunut 1900-luvulla ilmastonmuutoksen aiheuttaessa muutoksia suuren mittakaavan ja alueellisten ilmastoprosessien välisiin sisäisiin suhteisiin.

Eurasian large mammal dynamics in response to changing environments during the Late Neogene

Nisäkkäiden levinneisyyteen, niiden morfologisiin ja ekologisiin piirteisiin vaikuttavat ympäristön sekä lyhyet että pitkäkestoiset muutokset, etenkin ilmaston ja kasvillisuuden vaihtelut. Työssä tutkittiin nisäkkäiden sopeutumista ilmastonmuutoksiin Euraasiassa viimeisen 24 miljoonan vuoden aikana. Tutkimuksessa keskityttiin varsinkin viimeiseen kahteen miljoonaan vuoteen, jonka aikana ilmasto muuttui voimakkaasti ja ihmisen toiminta alkoi tulla merkittäväksi. Tämän takia on usein vaikea erottaa, kummasta em. seikasta jonkin nisäkäslajin sukupuutto tai häviäminen alueelta johtui. Aineistona käytettiin laajaa venäjänkielistä kirjallisuutta, josta löytyvät tiedot ovat kääntämättöminä jääneet aiemmin länsimaisen tutkimuksen ulkopuolelle. Työssä käytettiin myös NOW-tietokantaa, jossa on fossiilisten nisäkkäiden löytöpaikat sekä niiden iät.